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Abstract 



We study the heat equation with a random potential term. The 
potential is a one-sided stable noise, with positive jumps, which does 
not depend on time. To avoid singularities, we define the equation 
in terms of a construction similar to the Skorokhod integral or Wick 
product. We give a criterion for existence based on the dimension 
of the space variable, and the parameter p of the stable noise. Our 
arguments are different for p < 1 and p > 1 . 
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1 Introduction 



Our goal is to study the heat equation with stable Levy noise L(x) 
which depends only on space. 

Assumption 1 Throughout the paper, we assume that L has only 
positive jumps. 

Our noise is multiplicative, in the sense that it is multiplied by u in 
the equation. 

du . •, . .„ 

— = Au + uoL(x) (1.1) 
at 

u(0, x) = Uq(x). 

Here, L{x) is a one-sided stable noise of index p > 1 defined on x £ R d . 
The product u o L(x) is related to the Skorokhod integral or Wick 
product. We will give precise conditions later. Our main goal is to 
give conditions on p and the dimension d such that the solution exists 
and is unique. In fact, for p < 1, the stochastic integral is defined in 
a different way and can include more singular terms. In this case, we 
use the notation 

du 

— = A« + u*L(x) (1.2) 
u(0,x) = uq(x). 

There are many papers concerning the heat equation with a mul- 
tiplicative noise term (see |Par93j ) . In most cases, however, the noise 
depends on time, and is even independent from one time point to an- 
other. This is a very natural assumption from the point of view of Ito 
integration, since the relevant integrals often yield martingales. For 
example, the equation 

— = Au + f(u)W(t,x) (1.3) 
u(0,x) = uo(x) 

where W(t, x) is space-time white noise, is usually formulated in terms 
of an integral equation 

u(t,x) = / G(t,x - y)u (y)dy 

+ f f G(t-s,x- y)f(u(s, y))W{dyds) 

JO JR d 
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where G(t, x) is the heat kernel, and where the final integral can be 
defined in much the same way as the Ito integral. 

But if the noise does not depend on time, such an integral would 
be anticipating. Further evidence of the difficulty of noise depending 
only on space is given by the following equation: 

du 

— = Au + S(x)u (1.4) 
at 

u(0, x) = Uq(x) 

where 5(x) is the Dirac delta function. When u does not depend on 
t, this equation has been extensively studied by mathematical physi- 
cists. See the book AGHKH88 by Albeverio et al., for example. This 
equation (without t dependence) models the quantum mechanics of a 
particle under the influence of a point interaction. As usual, we can 
think of this equation in terms of the density of Brownian particles. 
The 5{x) term means that Brownian particles at would give birth 
to many new particles. But these new particles would themselves be 
close to 0, and so they might give birth to other particles, ad infini- 
tum. Because of this unstable behavior, when the dimension d is 2 or 
higher we need to take an approximate 5{x) function, multiply it by 
a very small term, and take the limit of the approximation. The case 
of x G R is better behaved. 

Yaozhong Hu has written two papers on time-independent noise, 
|Huf)lj and [Hu02 . He deals with the equation 

du . ■ , . 

— = Au + uoW(x) 
at 

u(0,x) = uo(x) 

where o denotes Skorokhod integration. One can also think of this 
product as the Wick product, see |Jan97j Chapter 3. W(x) is white 
noise in |Hu02| . but in HuOl it is colored Gaussian noise. By the 
Skorokhod integral, roughly speaking, we mean that in the relevant 
stochastic integrals, we should drop diagonal terms such as (W(x)) 2 . 
These terms correspond to repeated influence of W(x) for the same 
point x, so they are related to the singular behavior in the equation 

dm>. 

The goal of this article is to study for time-independent Levy 
noise L(x). The product o, like the Wick product, involves the deletion 
of diagonal terms in the appropriate integrals. But it is not defined in 
the same way as the Wick product for Brownian functionals, at least 
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with the usual definitions. Once again, the reason for using such a 
product is to avoid singular terms such as those coming from l|1.4|) . 
The product * used in (|l,2j) is similar, but we delete fewer terms. 

Our idea is that, once we understand the equation without the 
singular terms, we can try gingerly adding back the singular terms. 
But the understanding of the simpler situation must come first. 

In the time dependent situation, there have been several papers 
involving Levy noise, such as Kallianpur, Xiong, Hardy, and Rama- 
subramanian KXHR94 . In most cases, the noise has been additive. 
That is, it appears without any multiple of u or f(u). 

In the case p < 1, and for noise depending on both time and space, 
Mue98j dealt with the equation 

Flu 

^ = Au + u^L(t,x) (1.5) 

li(0, x) = Uq(x). 

If x G D for a smooth and bounded domain D C R d , and with some 
conditions on no, short-time existence was obtained in the case 

d< ^-P^ 
jp - (1 - p) 

For p > 1, Mytnik |Myt02| found that if 

2 

f>7 < "3 + 1, 
1 < p < min ^2, ^ + l^j 

then there exists a weak solution to (jl.5|) . 

Both of these papers were motivated by the superprocess with 
stable branching, see |Daw93j . We will not describe this superprocess, 
but it is believed that in some weak sense, it is associated with the 
equation Q1.5JI with 7 = 1/p. This motivation led us to study equations 
with one-sided noise. We leave possible generalizations to the reader. 

Our results are also divided into the two cases p < 1 and p > 1. 
Our main tools are multiple stochastic integrals, and this is the reason 
that we restrict ourselves to noise terms of the form f(u) o L(x) or 
f{u)-kL{x) with f(u) = u. Hu restricted himself to similar equations in 
the Gaussian case, for the same reason. Of course, multiple stochastic 
integrals have been used many times to study equations such as Ijl.HJl . 
For example, see Nualart and Zakai NZ89 and Nualart and Rozovskii 
NR97 . But in these articles, the noise always depends on time. 
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2 Theorems 



In this section we list our main results. They depend on some def- 
initions which will be explained later. We assume that uq(x) is a 
bounded function on D. Since our equation is linear, we may assume 
without loss of generality that ||t*o lloo < !• 
Here are our main theorems. 

Theorem 1 Suppose that < p < 1, and let uq(-) be a bounded 
function on D. Assume that either d < 4 or 

p<1 2 + dh> 

Then hl.ty) has a unique solution u(t,x). 

We are not sure that the conditions on p, d in the above theorem are 
optimal. 

Theorem 2 Suppose that p > 1 and uo(") * s a bounded function on 
D. Assume that 

2 

Then equation 11. 1)) has a unique solution u(t,x). Moreover, for all 
t > and x G D, u(t,x) G Lf oc (0), for all q £ [1,1 + (2/d)). 

Here L q loc (Q) denotes the space of all random variables / : ft — > R, 
such that, for all e > 0, there exists an event A C £1, of probability 
less than e, such that E[\f l^c \ q ] < oo. 

Remarks: 

1. In Theorem |21 in the Brownian case of p = 2, the critical dimen- 
sion is d = 2, which agrees with |Hu02j . 

2. For p < 1, if we use the product o instead of *, the proof of 
Theorem n is trivial. The reader can check that by removing the 
large atoms of the measure L(x), we can take q = 1 and use the 
proof of Theorem |2j 



3 Basic Definitions 

We will index sequences of real numbers, such as y,, by subscripts. A 
sequence of vectors in R n will be indexed by superscripts. Thus, x^ 
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is the i th vector in the sequence, and is the k th component of the 
jth vec t or _ 

If it is not stated otherwise, C in this paper will denote a constant 
which value my change from line to line. 

Throughout the paper, G(t, x) denotes the fundamental solution 
of the heat equation for x € R d , 



u(0,x) = 5(x). 

Explicitly, 

G(t ' X) = (4^ eXP ("¥)- 

Given a domain D C R d , we let G(t, x, y) be the fundamental solution 
of the heat equation in D, with Dirichlet boundary conditions. That 
is, G(t,x,y) satisfies 

An 

5 y (x) 

for x € 3D. 

We note the elementary inequality 

< G(t,x,y) < G(t,x-y). 

In other words, when G has two arguments, it is the heat kernel on 
R d , and when it has three arguments, it is the heat kernel on D. 
Finally, we write 

T n = T n (t) := {(si, ...,s n ):0<s 1 <...<s n <t}. 

4 The Noise 

Roughly speaking, our noise L(x) will be the analogue of dL(t)/dt 
where L(t) is a one-sided stable process of index p £ (0, 2), and 
t £ [0, oo). Of course, this derivative only exists in the generalized 
sense. For example, we could use the theory of Schwartz distributions. 

We will start with a review of the one-sided stable Levy processes, 
and refer the reader to Bertoin's book Ber96j for details. There is a 
basic difference between the cases p < 1 and p > 1. 



au 

~dt 

u(0, x) 

u(t, x) 
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For p < 1, we can construct L(t) via a Poisson process, as follows. 
Consider the following measure on (t, y) € [0, oo) 2 given by 

v{dydt) = c p y~ {p+1) dydt, 

where c p is a constant depending on p. This is the same constant 
that occurs in the Levy measure of the one-sided stable processes 
(see |Ber96j ). Let Y(-) be a Poisson random measure on [0, oo) 2 with 
intensity v. Then, the p-stable process L(t) is defined as 



p poo 

L t = / yY(dyds). 
J\o.t\ Jo 



'[o,t] Jo 

Note that if we let (U,yi) be the locations of the points (or atoms) of 
the Poisson measure Y, then 

In a similar way, we define the random measure L(A) : A C [0, oo) as 

l(a) = Y.y* 

ueA 



and observe that 



L t = f L(ds). 
Jo 



For p > 1, we must introduce an approximation and compensation 
procedure, as follows. Let 

Y n (dydt) = Y(dydt)l(y > 1/n). 
Note that if we define 

u n (dydt) := v{dydt)l(y > 1/n) 

then Y n becomes a Poisson measure with intensity u n . Note that 

v n ([0,t] x [0,oo)) =i/([0,t] x [l/n,oo)) = ^-n p t. 

P 

Next, let 

, . ft roo roc 

L t =\\ yY n (dyds) - 1 / ydu n {y) 
Jo Jo Jo 

t poo ~ 

\ yY n (dyds) - -^-nP-H. 
o Jo p-l 



There is a problem with p = 1. We leave the details to the reader, but 
roughly speaking, we would truncate the measure by removing large 
jumps. 

Note that L t is a martingale with respect to its own filtration. 
The p-stable process Lit) is the limit of the processes as n — > oo. 
For details of this limit, we refer the reader to |Ber96j . We could define 
the random measure L(dt) in a similar way. 

The notation L(t) stands for the the density of the measure L, even 
though this is only defined in the generalized sense. This is similar to 
the convention of writing B(t) for the derivative of Brownian motion. 

Our construction of L(A) : A C R d is similar. First consider the 
case p < 1. Fix a bounded open region D C R rf ; we will give some 
conditions on D later. For (x,y) G B C D x [0, oo), let Y(B) be 
the Poisson measure with intensity dxv{dy). Let (x^ l \yi) denote the 
locations of the points (atoms) of the Poisson measure Y. For A C D, 
we define 

L{A)= £ yi . 

xWeA 

For p > 1, we must approximate. Given n, let Y,v be defined as 
above, and let 

Y n {dxdy) = l(y > l/n)Y(dxdy) 
dxv{dy) = c p y~( p+1 ^dxdy 
dxv n (dy) = l(y > l/n)dxv(dy). 

As above, for A C D, we let 

L n {A) = ViHyi>Vn) 

= ViHvi > 1 / n ) 

where m(A) denotes the Lebesgue measure of A. Again, we leave the 
case p = 1 to the reader. To get L(dx) , we again take a limit as n — > oo. 
Since this procedure is not very different than for the one-sided stable 
processes with p£ [1, 2), we refer the reader to |Ber96j . 

For the future, we label the atoms of the measure L by (x^ l \yi) : 
ajW £ R d . Then, xW is the location of the atom, and yi is the mass 
of the measure L({x^'}). 



o J A 



ydxv n {dy) 



m(A) 



p — 1 



-TV 
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Fix K > 0. Since D is a bounded domain, there are only a finite 
number of atoms (x^\yi) with G D and ?/, > K. Let denote 
L with the preceding atoms deleted. That is, 

L K = L- ViS x (i) ■ 

Vi>K 

Assumption 2 In the succeeding sections, we will tacitly assume that 
L is actually Lk- Furthermore, let Ak denote the event that L = Lk, 
that is, there are no atoms larger than K . 

Note that 

lim P(A K )= lim P(L = L K ) = 1 

K^oo K— >oo 

From now on, we will replace L by Lk, and let K — ► oo at the end. 

5 Multiple Stochastic Integrals 

In the one-dimensional case, such integrals have been considered ear- 
lier, for example in RW86 . We suppose that D C H d is an open, 
bounded domain with smooth boundary. As described earlier, for 
p G (0,2), let L(A) be the one-sided stable random measure of index 
p, defined for Borel sets A C D. 

Assume that p > 1. Now we describe our multiple integrals with re- 
spect to the noise L, recalling that L(x)dx is just another notation for 
L(dx). Assume that we have a symmetric function f n (x^', . . . ,x^ n ') 
defined for x^ 1 ', . . . , x^ G D. Define D n C D n as follows. 

D n = {(xV),x&>,...,xV>)€D n : x? < x { ? < . . . < x^}. 

We regard x\ as the time variable, so that the following integral is 
defined in the Ito sense. 

Ufn)=n\ L f n (xW,...,x^)L(dx^)...L(dx^). 

Sufficient conditions for the existence of the integral I n (fn) are given 
in Lemma |21 

Since f n is symmetric, I n (fn) is invariant under permutations of 
the indices i = 1, . . . ,n of the x^ % >. For / = (/o, /i, . . .), where f n is a 
symmetric function of the variables x^\ . . . , x^ n \ let 

oo 
n=0 
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Here, /o is interpreted as a constant, and Io(fo) = /o- 

For future use, we define the symmetrization of a function as fol- 
lows. 

symCf)^ 1 ), . . . ,x(»>) = -i • • • 

n! ' 

where the sum is taken over all permutations a of the indices 1, . . . , n. 
If we wish to symmetrize over only some of the variables, we write 
those in a subscript. 

= ±J2f(x^,...,x^\zU,...,zW) 
n\ L — ' 

a 

Note that the integral /(/) has the following property. Because 
of the ordering we have imposed on the variables x®, it follows that 
1(f) depends linearly on the masses yi of the atoms of the measure L. 
That is, /(/) can have terms like yiy2, but no quadratic terms like y\ . 

We note in passing that our multiple stochastic integrals are re- 
lated to the Wick product for the compensated Poisson process Y; see 
Lytvynov Lyt03 , for example. 

For the case p < 1, we propose a different integral </(/), which 
allows for more singularities. Let X n (D) C D n be the set of coordi- 
nates (x^, . . . , x^) with no adjacent indices equal. In other words, 
x {i) x {i+l) for i = 1, . . . , n — 1. Let 

Jn(fn) = I f n (x^,...,X^)L(dx^)...L(dx^). 

JX n {D) 

Since L(dx^) . . . L(dx^) is a measure, the integral J n (fn) can be de- 
fined for each u, provided /„ is integrable with respect to this product 
measure. 

We define J(f) in a similar way: 

oo 

■/(/) = £ Ufn)- 
n=0 



6 A Framework for the Equation 

The goal of this section is to set up a framework for discussing the 
equation, so that existence and uniqueness can be rigorously discussed. 
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In the white noise case, there are some related concepts in NR97 and 
|MR98j . 

We define 

X n = {I n (fn) '■ fn is a symmetric function on D n }. (6.1) 

Definition 1 Let g(-, •) be a measurable function on D xQ, such that 
g(x, •) G I n for almost every x G D. That is, there exists a function 
f n on D n+1 such that 

g(x,-)=n\[ f n (x ( - 1 \...,x (n) ,x)L(dx ( - 1) )...L(dx^), a.e. x, 

and f n (-,x) is a symmetric function on D n for almost every x. 
Then we define 



g(x) o L(dx) 

(n + 1)! 1 S ym xm _ xin+1) f n (xW,. . . , ^,^ +1> ) 



D 

= [n 



L(dx^)...L(dx ( - n+1 ^). 

Set 

r oo i 
T = < y~] g n : g n G 2 n and the sum is converging in probability > . 
ln=o J 

Definition 2 Let g(-, •) be a measurable function on D x fi, such that 
g(x, •) G I for almost every x G D. That is, there exists a sequence 
{g n }n>o such that 

g n (x, •) £l n , n > 0, a.e. x £ D, 

and 

oo 

sOzv) = ^29n(x,-), a.e. x G D. 

n=0 

JTien we define 

/ g(x)oL(dx) = Y] / g n (x) o L(dx), 
Jd n=Q J D 

where the sum is converging in probability. 
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Next definition is analogous to the definition of Hu, from H uOlj . 

Definition 3 A measurable function u : R+ x D x Q i— > R is called 
a solution to the equation (1.1) if u(t,x) G I for every x G R rf and 
t G (0, T] and the following equation is satisfied: 



u(t,x) = / G(t,x,y)u (y)dy (6.2) 
Jd 

+ J D (J Q G (* ~ s ' x ' y) u ( s ' y) ds ^j o L ( d v)i 

Vx G R d , t G (0,T]. 
Our strategy is to expand u(t, x) in a recursively denned series: 

oo 

u(t,x) = ^2 u n (t,x), 

n=0 

where 

u (t,x) = / G(t,x,y)u (y)dy (6.3) 



u n+ i(i,x) = / / G(t - s,x,y)u n (s,y) o L(dy)ds. 
Jo Jd 

Note that Mo has two meanings. If uq(x) is a function of one variable, 
then it is the initial value for our SPDE If uo(t,x) depends on 

two variables, then it is the first term of the expansion for u(t,x). 
To deal with the case p < 1, we define the operation *. We define 

Jn = {Jn(fn) ■ fn is a function on £>"}. (6.4) 

Definition 4 Le£ g(-, •) be a measurable function on D xO, sitc/i i/iai 
<7(x, ■) £ J n for almost every x G -D. T/ioi is, i/iere exists a function 
f n on D n+1 such that 

g(x,-)= f f n (x^,...,x ( - n \x)L(dx^)...L(dx^), a.e. x, 
Jx n {D) 

and f n is a function on D n+l . 
Then we define 



D 



g{x) -k L(dx) 

f n (x {1 \. . . ,x( n+1 ))L(dx«) . . . L(dx( n+ V). 



>X n+1 {D) 

The other definitions for p < 1 are completely analogous to the 
case p > 1, with * instead of o. We leave this part to the reader. 
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7 Proof of the existence part of 
Theorem U (p < 1) 

7.1 Calculations for the heat equation 

As mentioned at the end of the Section |IJ we will assume that Ak 
occurs, so that L = Lk- We suppose that K = 2 N for some integer 
N. We will obtain conclusions that hold almost surely for each integer 
K = 2 N . Since P{Ak) —> 1 as K — > oo, our assertion will almost 
surely hold for L. 

We will use the notation of Mueller IMue98 . We call the atom at 
with mass ?/, a particle of type n if 

2 -("+l) < y. < 2 - n . 

Since there are no atoms of mass greater than K = 2 N , we need only 
consider particles of type n > —N. 

We will use the expansion defined in ()6.3|) . with * instead of o. 
Then uo(t, x) is the solution of the heat equation with initial function 
u (x). 

7.2 Distance between points 

We start with an elementary lemma about Poisson probabilities. 

Lemma 1 Let X be a Poisson random variable with parameter A. 
Then for any nonnegative integer uq we have 



P(X > n ) < 



n ! 



As a consequence 



A n ° 

P(X>n )<—e n °. 



n, 



no 




Proof. 

oo 

P(X > n ) = J2 P ( X 



n) 



n=no 
oo 

n=no 

oo 



A" 



A n A n-n n !(n-n )! _ A 



E 



"o! n t^ (n-no)! n! 
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But for all n > no, we have 

no!(n — no)! 1 



n! 



(-) 

\n J 



< 1. 



Thus, 



\no 00 \n— no 

P(X>n ) < ^ 



■e" A 



Since 



n o } - n t^ (" " wo) 
no! fr k\ 6 
n ! ' 



no Jin nn° 
= 1 + — + — + ... + — + .. . 
^ 1! ^ 2! ^ n ! ^ 



> 

we have 



n, 




n ! 



\«0 p no 

P ( . Y > „ o) < _ < r ._ 

I 

Suppose that D C B, where B is a cube in R d with side r > 1. 
Fix a,£ > 0. Consider the set of closed cubes of side length £ whose 
centers lie in the scaled lattice aZ d . Let Ce(B, a) be the subset of such 
cubes which intersect B. 

Suppose that x, y 6 D is a pair of points whose distance \x — y\ < a. 
We claim that there exists a cube in C3 a (B,a) which contains both 
x and y. Indeed, since the cubes C a (B,a) cover B, it follows that 
x € S a for some cube S a G C a (B,a). Let S3 be the cube with the 
same center as S a , but with side length 3a. It is easy to see that every 
point within distance a of S a lies in Ss a . Therefore, S^a is a cube from 
the set C3 a (B, a) which contains both x and y. 

Furthermore, let R C H d and let N(R, n) be the number of parti- 
cles of type n which lie in R. Let \R\ be the Lebesgue measure of R. 
Note that N(R, n) is a Poisson random variable with parameter 

\{R,n) = \R\ / u{dy) (7.1) 

J 2 -(n+l) 
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r 2 ~ n c 
= \R\ / ~^dy 

72-("+ 1 ) y p+ 

= C p \R\2 np . 

where C p = c p (2 p — l)/p. 

For a, M > 0, let A(n, m, M, S) be the event that within S Pi B, 
there is at least one particle of type n and at least M particles of type 
m, where S is a cube in C^ a {B, a). Let 

S = SDB 

and note that 

\S\ < {aArf 

where r is the side of B. We will estimate P(A(n, m, M, S)) for M = 
2 k , k > 0. There are two cases. 
Case 1: n ^ m. 

According to Lemma ^ 

P{A(n, m, M, S)) = P(N(S, n) > 1)P(N(§, m) > M) 

\{S,m) M e M 



< CeX(S,n 
= Ce\S\ M+1 



M M 

^np+mMp^M 



<ynp+mMp M 

< Ce.(aAr)W 2 ^ . 



Case 2: n = m. 



P(A(n,m,M,S)) = P(N(S, m) > M + 1) 

,A( g,n) M+1 e M 
M M 

nnp+mMp e M 

C\S\ M+1 



< c 



M M 

< C{ a,rf^^g^. 
Next, for some po > 0, let 

a = e - 1 l d fi i /{ Md )2~P n l( dM }2~ pm / d M 1 l d 2~ 5 ( n+m+k ^( dM ^ 
Note that 



dM _ c4 
l n,m,M ~ 



2'np+mMp e M 1 1 



M M 



2 



-S(n+m+k) 
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and therefore, for S G C 3an m M (£, a), 
(a n ^ 2k Ar)' d P(A(n,m,M,S)) 

^np+mMp gM 

< (a„, m , 2 * A r)- d (a„ imi2fc A r)^ M+1 ) ^ 

= ( a n,m,2 k A r ) 

< ^42-5("+m+fc) 

Let -A(-B) be the event that there exist n,m > — iV and M = 2 k with 
fc > 0, such that there exists a particle of type n in £?, and within 
distance a n ^ m ^M of this particle there are M particles of type m in B. 
Note that the number of cubes in Cs a (B, a) is the same as the number 
of cubes in C a (B, a). This number is less than or equal to a constant 
times (a A r)~ d . Therefore, we have the following. We write J2s f° r 
the sum over S G C^a k (B, a n m 2 fe). As before, S = S PI B. 

P(A{B)) 

oo oo 

< E E E p (^ m ' M ' s )) 



fc=0n,m=-jV 5 
oo oo 



= E E (v,2' Ar ) d E( a «rf Ar )^( i ( ti . m - M ' 5 )) 

k=0n,m=-N S 



oo oo 



< ^E E (a n , m , 2 ^r) d Y.^- 6in+m+k) 

k=0n,m=~N S 



oo oo 



= C£ E (a n , m , 2fc Ar) d 5 4 2-^ +m+fe )#5 
fc=0 n,m=—N 

oo oo 

= ^ E E <5 4 2~ <5 ( rH " rn+fc ) 




where C depends on N, but not on (5, as long as 5 < 1. Here, #5 
denotes the number of cubes S that we are considering. 
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Thus, we have shown: 



Lemma 2 There exists a constant C > such that the following 
holds. Fix a ball B C R rf of radius at least 1, let 5 € (0, 1), and let 

a = e - l / d ^/{ Md )2~P n l^ dM ^2~ pm ^ d M l ^ d 2~ 5 ^ n+m+k ^^ dM ^ 

Let A be the event that for some n, m, M , there exists a particle of 
type n in B with M particles of type m in B within distance a n m m . 
Then 

P{A) < CS. 



7.3 Estimation of the solution for the heat 
equation. 

Now we will use Lemma|2]to bound the solution u(t, x). Assume that 
event A does not occur, where this event was defined in Lemma |2j 

Taking x as fixed, define A\ to be the event that for some m, M, 
there exist M particles of type m within distance ao,m.M of x. The 
same arguments as before easily lead to the estimate 

P{Ay) < C5 

where C is the same constant as in Lemma |21 

For the rest of the section, fix n, and consider u n (t,x). Here n is 
no longer a subscript for a n ^ m ^i- We will assume that neither A\ nor 
A occur. Also, we write 

= X 

M = 2 k . 



Define 



M{n) = {(h, ...,i„)£N" : VI < k < n - 1, i k ^ i k+1 ) . 

Let v(t,-) be the solution of the heat equation on D with initial func- 
tion no, and recall that G(t, x, y) < G(t, x — y). Setting i n+ i = n + 1 
and m n+ i = 0, we get 

\u n (t,x)\l(A c nA c 1 ) (7.2) 

< E 

(ii,...,i„)6M(n) 
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n [G(s j+1 - Sj ,x^-x^) yi 



\v(si,x^)\ds 



< IKI 



E E 

...i„=Omi,...,m„=-iV 



n [g^ 

3=1 



J 1 rn,j+i,rnj,2 J ' 



ds. 



We assumed that ||tto||oo < 1- It is known that the solution of the 
heat equation is bounded when the initial function uq{x) is bounded, 
and \\v (t, •) ||oo < ||no||oo- When summing over all (ii, ... ,i n ) G M n we 
sum first over all possible types mi, . . . , m n of particles x^ n \ . . . , x^ l "\ 
For two fixed consecutive particles a;w and we let /cj be the 

smallest nonnegative integer k such that 



®"rn,j+i ,rrtj ,2 



ft < 



< a 



m.j +1 ,m.j, 



2 fe+l 



where m-j+i and m-j are the types of the particles x^i +1 ' and x^\ Ob- 



; riij+i aim 7ftj are me lypes ui me panicles x v J ri/ aim x v J '. wu- 
serve that such a fcj always exists. Because we are on the complement 
of both sets A and A\, there is no particle of type rrij at a distance 
less than a mj+1 ^ mj fi from x^ j+1 \ Thus, 



x 



> a m,j +1 ,m,j,2°- 



On the other hand, for fixed rrij+i and m.j, we have 



, lim a m m 2 k = 00. 



Thus for large enough fc, we have 



— a rnj +1 ,rrij,2 k - 



Therefore, 

G(a i+ i - s j ,a;fe+ 1 ) - x&>) < G(s i+ i - s^a^^*). 

Since x^ and are within distance a n fc,-+i for a fixed fe,- 

and a fixed x( l ' +1 ', there could be at most 2 k J +1 — 1 < 2 k i +l particles 
x^^ of type rrij within this distance from x^ j+1 \ This is the reason 
why the factor 2 k ^ +1 appears in the previous inequality. Since x^' is 
of type rrij, we have < 2~ mj and this is why 2 -m -> appears in (|7.2|) . 
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Estimating the term G in (|7.2|) . with M = 2 k , we get 
/ n 2 \ 

= Cr -d/2 c „„l U m 3 ,m j+1 ,M\ 



cxp 



Ar 



C r -^/2 exp ^ _ ^ r yl e -2/d s 8/(Md) 2 -2pm j+1 /(dM) 2 -2pm J /d 



^2/d2~2S(m j+l +mj+k)/(Md)j 



< Cr- d / 2 e* V ^-(Ar)- 1 e- 2 / d 5*l d 



22[k(M-6)-m J+1 (p+6)-m j (pM+6)]/(dM)j 



Therefore, 

G(r,a m . +1)m . )2fc ) < Cr^exp (.Cor-^^^ 1 -' 5 )-^^^/^^)]^) 
where Co depends on 5. Thus, 



u n (t,x) < 



E E 



n 

i=i 
• exp 



fci,...,fc„=0 m\,...,m n =—N 

C {s j+1 - Sj )- d/2 2 k i- m i 



(-C {s j+ i - s 3 -) _1 2 2 M 1 - 5 M m J + m J WM)( P +<5)]/^ 



(is. 



Rearranging the terms and using 2 < i an d since m ra +i = 0, 

we get 



u n (t,x) < 



n 

i=i 
■ exp 



C(s j+ i - s 



E E 

) m\,...,mn- 

d/22k j -(M/{M+l))(rn, j +m j+1 /M) 



Tn ki,...,k n =0m 1 ,...,m n =-N 



(7.3) 



(-Cq {s j+ i - Sj)" 1 2 2 [M 1 - 5 )-( m J+ m W M XP+' 5 )]A^ 



ds. 



We wish to replace the above sums over rrij by integrals. This would 
involve replacing each variable mj by a continuous variable Wj taking 
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values in [rrij, rrij + 1]. This would decrease the exponential, but would 
increase the first power of 2 by a multiplicative factor of 2 n . We would 
have 



\u n (t,x)\ < 

jTn k 1: ...,k„=0 

n 



V 



(7.4) 



n 

3=1 



C(s j+ i - Si )- d /2 2 fc j-( M /( M + 1 ))K+^+i/ A// ) 



• exp 



(-Co( 



Sj+l 



1 2^[kj(l-S)-(w J +w J+1 /M)(p+S)]/d 



dwds 



where, as with ds, we let dw = dw\ . . . dw n . Next, we wish to make 
the change of variables 

w j+1 



Vj = W-j + , r 

3 3 M 



J = l, 



, n. 



Now w n +i is fixed, and the reader can easily check that the Jacobian 
determinant for this change of variables is 1. Since M = 2 k > 1 and 
Wj > — N, we have that Vj > —2N. After making this change of 
variables, we can put the integral signs over v inside the product as 
follows. 



\u n (t,x)\ < 



E 



T - kl ,..., kri =o J - 2N 



2N 



(7.5) 



n 

• exp 



C( Sj+1 - Si )-d/2 2 fc;-(M/(M+l))^ 



(-Co( 



s i+l ~ s i 



1 2%j 0--fi)-Vj(p+5)]/d) 



dvds 



/ TT V / c( s , +1 - s ,)" d/2 2 fcj " (M/(M+1)) ^ 



j=l lkj=0 



• exp (-Co (s j+1 - sj)' 1 2 2 ^- s ^^/ d ) dv 3 



ds 



where once more, dv = dv\ . . . dv n . Let Tj^. denote the term inside 
the final sum in (|7.5|) , Dropping the subscripts on T,s,k and setting 
r = Sj + i — Sj, we can write 



Cr -d/2 2 k-(M/(M+l))v 



2N 



(7.6) 



• exp 



_ Cor -l 2 2[k(l-S)-v(p+S)]/d^ ^ 
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We will split the integral into two ranges for the variable v. 

For some values of v, the exponential in Q7.6JI is close to 1. This 
will happen if the terms following Cq in Q7.6|) are less than or equal to 
1. We must solve 

r -l 2 2[k(l-S)~v(p+5)}/d < L 
Simplifying the above expression, we get 

2 -2v(p+8)/d < r2 -2fe(l-<5)/d 



or 

M Md 

2-W+T V < r 2(M+l)(p+«) 2 (M+l)(p+i)' t _ (7_g) 



M(l-i) fc 



Let A = A(k,d,p) be the value of u for which equality is attained in 
(fTSjl . We will split T into two integrals, T = + T^ 2 ), where 



= / /(u)di; 
JA 

= / J(u)dv 



2N 

and is the integrand in (|7.6|) . 
Case 1: An estimate for T^ 1 ). 
By (|7.8j) . and using the fact that 

1 M 

- < < 1 

2 ~ M + l - 

we have 

^(1) = y°° (7 r -d/2 2 fc-(M/(M+l)> exp ^_ Cor -l 2 2[fc(l-*)-t;(p+i)]/^ dv 

/■oo 

< / (7r -d/2 2 fc-(M/(M+l)) V ^ 
J A 

< Cr -d/2 2 k-(M/(M+l))A 

d | Md .A M(l-<5) \ 

= Cr 2" l "2(M+l)(p+«) 2 V {M+l)(j>+«)/ 
j I d /./-i , M(l-i) S 

< (7 r a^4( P +«) 2 \ (m+ihp+,5) ; _ 

since < r < t. Here, C depends on i. 
If 

p < 1 
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and if 8 is sufficiently small, then (remembering that M = 2 A 

>H 1_ (M+l)(p+,5)J = T9 \ (2fc+l)(p+«5)y 



^ 2 fe l 1 -(M+l)( P +i)J = ^2 I (2>°+l)(p+6)) < 

k=0 k=0 



so that 



00 



E 2^<Cr"3 + W (7.9) 



Case 2: An estimate for T^ 2 ). 

Now we turn to the case v < A. Here, we assume that 

A > -2N 

or there will be no values of v for which — 2N < v < A. We will again 
drop the subscripts and set r = Sj + \ — Sj. Since equality is attained 
in ()7.7j) when v = A, we have 

r (2) 

c , r -d/2 2 fc-(M/(J 1 f+l))« exp /_ Cor -l 2 2[fc(l-<5)-.( P +5)]/^ du 

V / 

Cr-"/ 2 2 fc - exp (_c o2 2(A-*0(p+*)/<rt dv . 



-2AT 
A 



2N 

Changing variables to z = A — v, and since M/(M + 1) < 1, we see 
that 

r (2) = J^ 2N Cr -d/2 2 k-(M/(M+l))(A-z) exp (_ Co2 2z(p + 6)/d^ ^ 

< (7 r -d/2 2 fc-(M/(M + l))A y°° 2 z exp ^_ Co2 2 Z ( P+ 5)/^ ^ 

< Cr -d/2 2 k-(M/(M+l))A 

< (7 r 2 + 4(p+6) 2 V 1 (M+l)(p+«) J _ 

Once again, if 

p < 1 

and 5 is small, then 

00 ,/, M(l-<5) \ 00 ,Y X 2 fc (l-5) \ 

fc=0 fc=0 
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so that 

T (2) < Cr"^ +3 (^). (7.10) 
Combining (f7T9|) and (|7~T0|) . if 

p < 1 

and 5 is sufficiently small, then 

T = + < Cr~^ +T <^> 

and 



u n (t,x) <C H (s i+1 - si)-~-~—^ (1*. 



Now, a calculation of Hu ( |Hu01| . proof of Theorem 4.1, equation 
(4.10)) gives, for > 0, 

n /n(l-)3)rn/i _ a\ 

T.n^ = r ( i + ,(L ffl T <™> 

We give a quick derivation of Hu's estimate in the appendix, since he 
merely cites the result. 
In (17.111). substitute 



2 4(p + 5) 

If f3 < 1, then T(l — /3) is defined. With C n denoting a term of order 
C n for some C > 0, we get 

T(l + n(l-/3)) 

since only the denominator has order greater than C n . Thus, if (3 < 1, 
and if A and A\ do not occur, 

oo 

u(t,x) = ^ u n (t,x) < CO. 

n=0 

Solving for (3 < 1, we obtain the requirement that 

d d 

2 ~ 4(p + 5) < 
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For small enough values of 8, this is satisfied if either d < 4 or 

This completes the proof of Theorem^ for Lk- As mentioned earlier, 
since P{Ak) = P(L = Lk) — > 1 as K — > oo, it follows that Theorem 
^ holds almost surely. 



8 Proof of the existence part of 
Theorem El (p > 1) 

For p > 1, we also wish to replace L by L/<. But in this case our 
stochastic integrals will no longer be martingales unless we delete the 
compensator for L — Lk- We define this compensator as follows. 

C(dtdx) = dtdx J c p y- (p+1) dy = c' p K- p dtdx. 

Then, on the event Ak, we can replace (jl.lj) by 
du 

— = Au + c' p K~ p u Lk(x) 
u(0,x) = uo(x). 

If we let 

w(t,x) = exp ^— c p K~ p tj u(t,x) 

then w(t, x) solves 

du ■ , . , 

— = Au + uoL K (x) (8.1) 

u(0, x) = uq(x). 

Therefore, any conclusions which hold almost surely for w(t, x), for all 
K, must hold for solutions to (jl.lj) . 

For ease of notation, we write u(t,x) instead of w(t,x) and L in- 
stead of Lk- With this notation, we can prove the following lemma. 



Lemma 3 For any q G (p, 2] there exists a positive constant C q such 
that for any measurable symmetric function f : D n -*Rk have: 

f(x^\. . . , x^)L{dx^) . . . L{dx^) q 



E 



IV. 



< 



C"(n!) 



D n 



\f(x^,...,x^)\ q dx^ ...dx (n \ 
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Note that the above lemma gives sufficient conditions for the existence 
of the multiple stochastic integral I n (f n ). 

Proof. We would like to use the Burkholder-Davis-Gundy inequal- 
ity, and for this reason we need to expand on our formulation of the 
multiple stochastic integral given in Section (JHJ) , so that we have a mar- 
tingale. We start with the a-algebra. Recall that we have chosen the 
first coordinate of x as our time variable. In particular, we can define 
a filtration. Let m(D),M(D) be the minimum and maximum values 
of xx, the first coordinate of x, for x £ D. For z E [m(D), M(D)\, we 
define the cr-algebra 7i z as follows. 



n z = ct{l(A) :Ac{xeD: Xl <z} 

Next, let 



R(z,D) = {x € D : x\ < z}. 

Let h(x) be a predictable function with respect to the TC Z filtration. 
Let 



Y z = / h(x)L(dx). 

JR(z,D) 

Then (Y Z ,H z ) is a martingale, for appropriate conditions on h, and 
its quadratic variation is 

(Y) z = J2 h\x^)yl 
For ease of notation, we drop the subscript z when z = M(D), so that 

(Y) = (Y)m(d)- 

The Burkholder-Davis-Gundy inequality now implies that for q > 1, 
there is a constant C q depending only on q, such that 

E[\Y\ q ) < CgE \(Y) q ^ 

Recall that for r < 1 and a sequence of nonnegative numbers dj, the 
following elementary inequality holds. 

(oo \ r oo 
j=l / j=l 
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Setting r = q/2, and using our formula for (Y), we get 

' oo \ <?/2 



E[\Y\ 



< CE 



< CE 



\j=i 

oo 



Eh 2 (- (4) ) 

i=l 

C [ E\h(x)\ q dx 

JD 



1)1 



q/2 



(8.3) 



provided q > p. The constant C depends on p, q, and K. Applying 
(|8.3f) repeatedly, we get 



E 



D, 

< c n 



f(x ( - 1 \...,x ( - n) )L(dx W )...L(dx ( - n) ) 
f{x (1 \...,x (n) ) 9 dx W ...dx {n) 



(8.4) 



Dn 



Note, that the repeated application of inequality (|8,3|) was possible 
due to the order Xi < Xi < . . . < x± that we imposed on D n . To 
see this, let us remember that the function h involved in (|8,3|) needed 
to be predictable. Thus, 



E 



D, 
E 



f(x^\...,x^)L(dx^)...L(dx^) 



D 



h(xW)L(dxW) 



where 

h(x (n) ) 

and 

#n-i (2/) 



/(x«, x {n ^\x {n) )L{dx W ) . . . Lidx^-V) 



X 



(1) 



, X 



(n- 



1} )e^ 



_l : x 



(n-1) 



<y\ 



(8.5) 



where y £ R. The function h(x^) is predictable since x^ 1 < a?]** , 

we can apply the inequality (|8.3|) 

and obtain 



E 



L f{x {1 \. • • ,x^)L{dx<V) . . . L(dx^) 
Jd„ 



< C E 

JD 
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and so on. Using l|8.3|) . and the fact that / is symmetric, we get: 

11 



E 



n\ I f{x w ,...,x w )L{dx w )...L{dx w ) 



(n\) q E 



f(x^,...,x^)L(dx^)...L(dx^) 



< (n\) q C n I fix 



{n\) q C n 



1 



dx {l) . . . dx {n) 



(n)> 



C^n!)^ 1 / f(x- r 



(1) Jn)^ 



' dx (1) . . . dx {n) 
...dx {n) . 



Before proving Theorem [2 we present two simple lemmas. 

Lemma 4 Let n be a positive integer and q a real number, such 
that q > 1. Let f : D n — > R be a function in ~L q (D n ,dx), where 
dx is Lebesgue measure. If f is the symmetrization of f , then f € 
L q (D n ,dx) and 

\\f\U < 11/11. (8-6) 

where \\ ■ \\ q denotes the norm of the space h q (D n ,dx). 

Proof. For any permutation a of the set {1,2, ... ,n}, we can 
consider the function f a : D n — > R, defined by 



Ux^,...,x^)=f(x^,...,x^>) 

, y 



Making the change of variable g/ 1 ) = x^ 1 )) 
can see that 



(n) _ x {<r{n)) 



we 



D" 



\f a (x^,...,x^)\ q dx ( - 1 K..dx^ 



\f(x^,...,x^)\ q dx^ ...dx( n \ 



Thus for all permutations a, we have 

II fa \q = || / \\q ■ 

Since the symmetrization / of / is defined as 

n! ^ 



28 



we can apply the triangle inequality to get 



fo- \\q — II . ^fcr 



< 



\ ii f° 1 1 9 



1 



Ell/ 
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In the proof of Lemma Q we used the inequality e n > n n /n\, for 
all positive integers n. We would like to prove this inequality for all 
positive real numbers x and therefore we have to replace n! by 
Thus we will prove the following lemma. 

Lemma 5 For any positive real number x, we have 

+ > x x e~ x . (8.7) 

Proof. If x > 0, then for all t > x, t x > x x . Thus we have: 

coo 



poo 

T(x + 1) = / t x e- l dt 
Jo 

poo 

> / t x e- l dt 



oo 



> / x x e~ t dt 

J X 

poo 

= x x / e-*dt 



We can prove now Theorem |5J According to inequalities (|8.2|) and 
(|8.6j) . we have: 



E [\u n (t,x)\ 



91 



< C^niy- 1 / sym x(1) s( „)|t£o(* (1) )| 
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dx 



dx, 



„ n 

■ / ]jG(s k+l -s k ,x^ -x^)ds 
jTn k=i 

< aC n {n^- 1 f f f[G(s k+1 -s k ,x( k +V -xW)ds 

J D n JT n k= ^ 

where a =|| uq ||oo and 

T n = {(si, . . . , s n ) G R n : < si < s 2 < . . . < s n < t}. 



We may assume that || uq ||oo — 1. Let V(T n ) = t n /n\ be the volume 
of T n . Using Jensen's inequality, we obtain 



. n 

/ Y[G(s k+1 -s k ,x^-x^)ds 
jTn fc=i 

„ n 

n fc=l 

„ n 

< ^(TnJ'yCTn)- 1 / J] G"(s k+1 - s k ,x^ - xW)ds 
jTn k=l 

j.n \ q—1 f n 

' I l[G«(s k+1 -s k ,x( k+ V-xW)ds. 

T " k=l 



1.9) 



Recall that 

G q (r,x)dx = Cr-ti-W 2 J j~ d l 2 exp 



R<< 



R<< 



qx 
4r 



dx 



(8.10) 



Combining (|8.8|) and (|8.9|) . we see that (n!) 9 1 is canceled out by 
(1/n!)" 7 " 1 , and by using Fubini's theorem and Q8.1UJI repeatedly (i.e. 
integrated first over xW, then x^ 2 \ and so on), we obtain: 



E[\u n (t,x)\ c ' 



.11) 



< (ct^Y f [ f[ Gi{s k+1 - s k , x^ - zM)dx 

< (ct^Y [ [ ft G 9 (s fc+1 - Sfcj - x«)dx 

„ n 

= {ct^yj iKs^-s^-^ds. 



ds 



ds 



fc=i 
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Now, an estimate of Hu (|7.11|) gives 



k=l 



t n(l-{d/2){q-l)) T n^ _ |(g_ 1)) 

r(i + n(i-£( g -i; 



If d < 2/{q - 1), then T(l - (d/2)(q - 1)) is defined. For d < 
2/(q - 1), let a = 1 - {d/2){q - 1) > 0. Then from (I8TTT) we conclude 
that 



E[\u n (t,x)\ q ] < 



.12) 



T(l + na.y 

for some constant p depending on t, q, p, and K. Using the well- 
known inequality (b + c) q < 2 q ~ 1 (b q + c q ), for all q > 1 and 6 and 
c non-negative numbers, an easy induction argument shows that for 
a n > 0, and c = 2* -1 , 

/ oo \ 9 oo 



\n=0 



n=0 



Applying this inequality to u(t,x), we find 



E[\u(t,x)f 



E 



< E 



\ u n(t,x)\ 
\n=0 ) . 

oo 

J2c n+1 \u n (t,x)\ q 
Y,c n+1 E[\u n (t,x)\ q ] 

(c P r. 



$.13) 



n=0 

oo 



- C 5 r(l + na) 
Using now (|8.7|) we obtain: 

E[\u(t,x)\«) < cJ2 



n=0 



(na) T 



$.14) 



Thus, we can see from (j8.14|) that if <i < 2/{q — 1) (or equivalently 
q < 1 + (2/d)) and p < q, then £ [\u(t, x)\ q ] < oo. (The condition 
p < q was used in the Lemma 0) For such a q to exist we need that 
p < 1 + 2/d. 
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Let II be the random Poisson set {(x^ ,yi)}i>±. For any positive 
integer k we define the event: 

A k = {lo G SI : V (x,y) € U(uj),y < k}. 

That means, is the event that all atoms have masses less than or 
equal to k. 

Thus if p < 1 + (2/d), then for all g G (p, 1 + (2/d)) and for all 
k>l, 

E[\u(t,x)\l q Ah ]<oo. 

Since, the expectation is the integration with respect to a probability 
(finite) measure, it follows, that for any q' and q, such that 1 < q r < 
p < q < 1 + (2/d), we have: 

(i?[|n(t,x)U fc r']) W < (^[K(i,x)U ft r])^ < oo. 

Thus the solution u(t,x) exists, for all t > and all x £ D, and 
belongs to the space L^ oc (r2), for all 1 < g < 1 + (2/d). 

This completes the proof of Theorem |^1 if u(t,x) means w(t,x) 
and L means Lk- As mentioned at the beginning of Section El this 
proof carries over to u(t,x) and L, since P(Ak) = P{L = Lk) — > 1 
as if — > oo. 

9 Uniqueness 

9.1 The Main Lemma 

We consider both cases p < 1 and p > 1. Suppose u(t, x) is a solution 
to (jl.lj) in case p > 1 or H1.2|) in case p < 1. Then we have 

oo 

u(t,x) = ^2g n (t,x) (9.1) 

n=0 

where g n € in case p < 1 or g n G X n in case p > 1. 

Lemma 6 Letu(t,x) be as in \9.1\) . If the g n are uniquely determined 
from u, then uniqueness holds for hl.l\) or \1. 6 J\) . 
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Proof. Both cases have the same proof, so we will only deal with 
the case p > 1. Suppose there are two solutions Ui(t,x) for i = 1,2, 
where 

oo 

Ui(t,x) = 9n\t,x) 
n=0 

and gn\t,x) € 1 n . Let 

oo 

Y(t,x) = ui(t,x) -u 2 {t,x) = ^2g n (t,x) 

n=0 

where g n (t,x) = g2 (t,x) — gn\t,x). Then Y(t,x) is also a solution 
to with initial condition Y(0, x) = 0, so by ()6.2|) . we have 

oo 

^2g n (t,x) = Y(t,x) 

G(t - s, x, y)Y(s, y)ds^j o L(dy) 

^ Id Uo G(yt ~ s ' x ' y^ 911 ^ 7 y ^ ds ) ° L ^ 



n=Q 



D \J0 



n=0 ' 

oo 



Vx G R d , i € (0,T]. 



n=l 



where 



g n (t,x) = J d (^J q G{t - s,x,y)g n -i(s,y)ds ) j o L(dy) (9.2) 

Since we are assuming that Y(t,x) uniquely determines both g n (t,x) 
and g n (t,x), it follows that go(t,x) = and 

g n (t,x) = g n (t,x) (9.3) 

for any positive value of n. Thus, (|9.2|) implies that gi(t,x) = 0, 
and hence by H9.3|) g\(t,x) = 0. An induction argument now yields 
g n (t, x) = for all values of n. I 
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9.2 Proof of Uniqueness for p < 1 

For simplicity in this proof we will assume from the beginning that L 
is a stable noise (not Lk) without truncation of large atoms. 

We begin by ordering the atoms of L according to size, so that 

yi>V2>---- 

For z £ (1)2), let Ri(z)L be the transformation of the noise L 
which replaces yi with zyi. 

Lemma 7 Fix z 6 (1,2), i > 1. Let 

M = Ri{z)L 

Also, let Pm,Pl be the probability measures on the canonical proba- 
bility space of atomic measures, generated by M and L, respectively. 
Then 

Pm «Pl- 
Proof. For any k > 1, n > define an event 

Ak, n = |there are exactly n atoms of L 

with masses greater than or equal to 1/k} 

Note that for any k > 1 

oo 

(J A k,n = « 
n=0 



It is obvious that 



i 



lim P L |J A ktl = 0. (9.4) 
k ^°° \i=o J 

Fix 5 > arbitrary small. By (|9.4|) we can fix k sufficiently large such 
that 

PL[\jA k A<5. (9.5) 
\i=o ) 

Define 

/ i \ c oo 

A k)i = |J A k>l = \J A k>l . (9.6) 

\l=0 / l=i+l 
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Now will show that 

P m «Pl 0114,. (9.7) 
Let us decompose L as follows: 

L = L^ k + L 2 ' fc 

where L\ includes atoms of size y > 1/k and L2 includes atoms of size 
y < 1/k. Note that L l,k ,L 2,k are independent, and 

Ri{z)L = Ri(z)L^ k + L 2 ' fc onA k>i . 

Hence, if we define 

Ml = Ri(z)L^ k 

then, in order to get (|9.7|) . it suffices to show that 

P Ml « P L i,k on A k>l , VZ > i + 1 (9.8) 

where PmxiPl 1 * are the measures induced by Mi,L l,k respectively. 
Fix an arbitrary I > i + 1. On , Mi and L 1, have Z atoms at the 
same locations. Therefore, to get the desired absolute continuity we 
just need to show the absolute continuity of the laws of the correspond- 
ing atom sizes. To be more precise, from now on we assume that A k i 
occurs. Let Y-y , . . . , Yj (resp. , . . . , Yj M ) be the unordered sizes 
of the atoms of L 1,k (resp. Mi). {Y\ Y t L ) and (Y^ Y t M ) 
are continuous /-dimensional random variables whose laws are sup- 
ported on (l/k,oo) 1 . Moreover, the probability density function of 
(Yi , . . . , Yj ) is positive at any point of (1/k, oo) 1 . This immediately 
implies that conditioned on A k i the probability law of (Yj , - • • , Y t M ) 
is absolutely continuous with respect to the law of (Yj , . . . , Y, L ), and 
since I > i + 1 was arbitrary, ()9.8|) and ()9.7|) follow. 

The end of the proof of the lemma is trivial. Let B be any event 
such that 

Pl OB) = 0. 



Then 



P M (B) = J2 P M(Bf]A ktl )+ J2 PM(Bf]A k)l ) 

1=0 l=i+l 



J2 P M(Bf]A 



fcjj 



1=0 



< P M (\jA k>l ), (9.9) 
1=0 
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where the second equality follows from Q9.7|) . Now it is easy to check 
the following equality of events: 

| there are no more than i atoms of L with masses y > l/k\ 
= | there are no more than i atoms of M with masses y > l//c| 

Therefore, we get 

Pm (\jA k A = P L [\jA k \ . 

\l=0 J \l=Q J 

Thus, and (JHHJ) imply that 

Pm(B) < 5. 

Since 5 > was arbitrary, we get Pm{P>) = and the result follows. 
I 

By Lemma |H1 the following lemma gives the required uniqueness 
for the case p < 1. 

Lemma 8 Let p < 1. 7n HOl^ ^ e functions g n are uniquely deter- 
mined. 

Proof We argue by contradiction. Suppose Lemma El were false. 
Then by subtracting the two series, we form a series 

oo 

9 = E 9n = 

n=0 

where g n G J n , and at least one of the g n is not identically 0. In fact, 
we may assume that 

9n = ^ t Vi\ ■ ■ ■ Ui n n n i ■ ■ ■ 3 x i n ) 

(ii,...,i n )eJn 

with 

00 

E E • • -yiJhniXh,- ■ ■ ,%i n )\ < 00 (9-10) 

n=0 (ii,...,i n )£j n 

with probability 1. Actually, each solution is a difference of two series 
for which the functions h n are nonnegative. 
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Recall that {yi > 1/2 > ■ ■ ■} are the ordered sizes of the atoms of 
our noise L. Fix an arbitrary n > 1. Our goal is to show that, almost 
surely, 

Vi ■ ■ ■ y n K(xi, ■■■ ,x n ) =0 

from which it follows that h n = almost everywhere. By Lemma[3 for 
z; L £ (1, 2), the probability induced by Ri(zi)L is absolutely continuous 
with respect to the probability induced by L. Furthermore, Ri(zi) 
extends to an operator on J , as follows. We set 

00 00 

Ri(zi) ^29n = Yl R ii z i)9n, 
n=0 n=0 

where we define Ri(zi)g n to be a multiple stochastic integral with L 
replaced by Ri(zi)L. Therefore, if Q = J2^=o9n = almost surely, 
then Ri(zi) Y^=$9n = E^o R i{ z i)9n = almost surely as well. Fur- 
thermore, if Zi € (1,2), then 

Ri(zi) . . . R m {z m )g 

is an analytic function in z\, . . . , z m . Let A ntm Q be the terms in the 
expansion of this analytic function which contain Z\ . . . z n but no other 
Zi. Note that each A ntTn Q equals almost surely, and that 

lim A nm g = y 1 . . .y n h n (xi, . . . ,x n ) = 

m— >oo 

almost surely. In fact, the convergence to follows from ()9.1U|) . since 
the operator A„ )m removes more and more terms as m increases. This 
proves Lemma El I 

9.3 Proof of Uniqueness for p > 1 

By Lemma H3 the following lemma gives the required uniqueness for 
the case p > 1. 

Lemma 9 Let p > 1. In H9.1\) . the functions g n are uniquely deter- 
mined. 

Proof. Fix (t,x) and suppose that u\(t, x), U2(t, x) are two solu- 
tions. Subtracting, we obtain 

Y(t,x) = U\(t, x) — U2(t, x) 
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oo oo 

= £<#>(t,s)-2>«(t,*) 

ra=0 n=0 
oo 

n=0 

where gn (t,x) corresponds to the solution Ui(t,x). Also, each g n can 
be expressed as a stochastic integral 

g n (t,x)= L f n (t,x;x( 1 \...,xW)L(dxW)...L(dxW). 

Adding these terms, and using the notation from the beginning of 
Section |HJ we find that 

Y(t,x) = fo(t,x)+ [ h{t,x;y)L(dy) 

JR(M(D),D) 

where 
fx(t,x;y) = fi{t,x;y) 

oo - 

+ E /- /»(*> ^ ^ • • • > x [n - l \y)L{dx^) . . . L(dx^) 

and R n -i(y) was defined in (|8.5j) . 
Since x) = and 

z -> / f 1 (t,x;y)L(dy) 

JR(z,D) 

is a local martingale, it follows that fo(t,x) = and that the 
integrand fi(t, x;y) = for almost every value of y. But fi(t, x;y)=0 
is itself a sum of a deterministic function and a stochastic integral. 
Therefore we may use the same argument to show that fi(t, x;y) = 
for almost every value of y. Continuing, we can use induction to show 
that for each value of n, f n (t, x; y\, . . . , y n ) = for almost every value 
of yt,...,y n . I 



A A Derivation of Hu's Estimate 

In the proof of Theorem 4.1, equation (4.10) of Hu02 , Hu gives the 
following estimate (in his equation, a = 1 — d/4, and we use (3 = d/A). 

„ n 

Knit) := / H(s k+ i -s k )- p ds. 
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Then 

(Hu uses (3 = 7/4. Also, Hu has 2 instead of 1 in the above formula. 
His answer is equivalent if we ignore multiplicative factors of n.) 
Making the change of variables s = ar, ds = adr, we find 

„ n 



/ T ^n( 

= a - nl3 a n n n (t) 

= a n ^H n (t). 



Then we get the recurrence 

ft 
10 



H n+1 {t) = f (t - s)- p H n {s)ds 
Jo 

= [\t - S)~P S^-^Hn^ds 

Jo 

= H n (l) l\t- S )~^S n ^^ds 

Jo 

= t- n ^-^n n {t) f t {t-s)- f3 s n 

Jo 



Let 



a = (3 

b = n(l-/3). 

Recall that the Euler beta function is defined as follows, for a, b > 0. 

B(a,b) := [* (1 - xf^x^dx. 
Jo 

It is a fundamental identity that 

Changing variables, we get 

Jo [ } T(2-a + b) 



39 



Therefore, 



n n+1 (t) = t-^-^Hn^t 1 -^ 1 -^ 



r(i-/?)r(i + n(i-/?)) 

r(2 - P + n(l - 0)) 

r(i-/3)r(i + n(i-/?)) 

r(l + (n + l)(l-/3)) 



The product telescopes. Since 7ii(t) = t 1 we get (IA~T1) . 
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